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SUMMARY 


An investigation has been conducted to determine minimum "maneuver costs" 
for attitude reorientation of spacecraft of all possible inertial distribution 
over a wide range of maneuver angles by use of a two-impulse coning method of 
reorientation. Maneuver cost is proportional to the product of fuel consumed 
(total impulse) and time expended during a maneuver. Assumptions included 
external impulsive control torques, rigid-body spacecraft rest-to-rest maneu- 
vers, and no disturbance torques. 

Results are presented in terms of average cost and standard deviation for 
various maneuver ranges. Costs of individual reorientations can be calculated 
with the computer program included. 


INTRODUCTION 

Two-impulse attitude reorientation is a method of reorientation in which 
torque impulses are used to initiate and (later) to terminate free precessional 
and nutational motion of a body. Between torque impulses, the spinning body 
coasts undisturbed through some desired attitude change. The method has poten- 
tial application to maneuvering spacecraft capable of producing directed exter- 
nal torque impulses. 

Attitude reorientation by means of torque impulses has been investigated 
by several authors (refs. 1 to 4) . Reference 1 presents a two-impulse reorien- 
tation scheme for large-angle spin-axis reorientations of spinning axisymmetric 
bodies. Spin-axis or pointing-axis reorientations are defined as general 
reorientations with final attitude being arbitrary about the spin or pointing 
axis. Reference 2 extends the spin-axis reorientations of reference 1 to asym- 
metric bodies. However, precession angle is restricted to 180° and spin-axis 
inertia must be larger than inertias about transverse axes. Reference 3 intro- 
duces a statistical "average cost" for general reorientations of nonspinning 
bodies over a specified maneuver-angle range and compares two— impulse reorienta- 
tion costs with costs of other maneuver schemes. Results are limited, however, 
to elongated axisymmetric bodies. Reference 4 extends previous results to 
include two-impulse general reorientation costs for nonspinning axisymmetric 
foreshortened bodies. Also, two- impulse reorientation costs are obtained for 
pointing-axis reorientations of asymmetric nonspinning bodies of all possible 
inertial distribution. Finally, reference 4 presents three-impulse general 
reorientation costs for asymmetric nonspinning bodies. These maneuvers con- 
sisted of a basic precessional-nutational motion followed by an impulsive spin 
maneuver. This motion combination was required to satisfy an arbitrary con- 
straint that initial and final angles between the spin axis and momentum vector 
be equal. 

The present paper also presents general reorientation costs for nonspinning 
asymmetric bodies. However, the aforementioned motion constraint was eliminated 



through the use of elliptic functions and a computational scheme involving 
nested Newton-Raphson iterative loops. As a result, general two- impulse solu- 
tions involving precessional-nutational motion were found with significantly 
reduced costs. 

In the present investigation, reaction thrusters are assumed to initiate 
and terminate all maneuvers with impulsive torques of negligible duration. 

Other assumptions include rigid-body spacecraft, lack of disturbance torques, 
and zero initial and final attitude rates. 

The product of total maneuver impulse (proportional to fuel consumed) and 
total maneuver time is determined for the various solution paths available for 
each maneuver. This product is then nondimensionalized by spacecraft inertia 
about the intermediate axis and is referred to as the cost function for that 
maneuver. Since total impulse and maneuver time are inversely proportional for 
impulsive maneuvers, results are independent of each of these quantities. The 
desired or optimal solution is defined as the solution associated with the 
smallest cost function for each maneuver. 

Since it was assumed that all reorientations within a given range of maneu- 
ver angle are equally probable, optimal costs for a large number of statistically 
representative reorientations within a given maneuver range were averaged to 
represent reorientation costs for that maneuver range. Maneuver ranges of 
tt/360, tt/4, tt/2, 3tt/ 4, and tt radians were used. These results are intended 

primarily for reorientation cost comparisons of the two- impulse coning method 
with other reorientation methods. However, such comparisons are not included 
in the present paper. 

Although costs for individual reorientations are not presented, such costs 
can be determined, if necessary, through the use of the computer program pre- 
sented in appendix A. 


SYMBOLS 


cn (u,m) 

dn(u,m) ; Jacobian elliptic functions (see eqs. (12)) 
sn (u,m) J 

F nondimensional asymmetry factor defined by equation (9b) 

H total angular momentum 


Hx/Hy,H2 

^xy 

Ix/^y^^z 


angular momentum components along x-,y-, z-axes, respectively 

angular momentum in transverse plane, + Hy^ 

principal body inertias about x-,y-, z-axes, respectively; ly is 
intermediate manent of inertia 


m 


parameter of Jacobian elliptic functions defined in equation (14e) 



P,QfR 


motion constants (see eqs. (14a) to (14c)} 


Px 

Q2 ' 

Q2,ArQ2/ 

Q2,sw 

R) fR2»R3 

^2, max 

P2C 

SDA 

s 

T 

t 

u 

X,Y,2 

x,y,z 

I 

I 

X 

$ 


inertial factor (see eq. (27)) 
cost function (see eq. (26)) 

axial and transverse components of Q 2 , respectively (see eqs. (25)) 

product of total impulse and total time divided by ly for coning 
maneuver where thrusters, on spin axis can be swiveled 

sequential Euler angle rotations about z-,y-,z-axes, respectively, 
relating body axes before and after reorientation (see fig. 4) ; 

R 2 is also referred to as maneuver angle 

maneuver range, rad (0 5 R 2 ^ P2,max^ 

maneuver cuigle transformed for x-axis maneuver 

ratio of standard deviation to average cost 

1 

= -Oo + R2 + 0f) 

kinetic energy (see eqs. (15)) 
time, sec 

argument of elliptic function defined in equation (13) 
inertial coordinate axis system 

body principal coordinate axis system; unless otherwise specified, 
z is spin axis 

Xf-axis translated fran mass center along Zf-axis 

XQ-axis translated from mass center along ZQ-axis 

constant of linearity (see eqs. (13) and (14d)) 

initial base angle of spherical triangle, used to start iterative 
solutions, IT - R-j - (|)q 

rotation sequence of Euler angles about z-,y-,z-axes, respectively, 
relating inertial and body coordinate axis systems 

precession rate, nutation rate, and spin rate, respectively; referred 
to as Euler rates 

body rates about x-,y-, z-axes, respectively 
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Subscripts: 

av average value (with respect to time) during maneuver 


f 

value 

at finish of coning maneuver 

max 

maximum value 

new 

value 

from present iteration 

o 

value 

at start of coning maneuver 

old 

value 

from previous iteration 

start 

value 

to start iterative solution 


ANALYSIS 

Equations of Motion 

Two reference frames used in the analysis of the impulse coning method of 
attitude reorientation are illustrated in figure 1 . The body principal axis 
system x,y,z is related to an inertial axis syst«n X,Y,Z by an Euler rota- 
tion sequence about body axes z,y,z, respectively. Note that 4* and 

Z are chosen to be in the same direction as the total angular momentum H. 
Angular manentum H is established for each reorientation maneuver by the ini- 
tial torque impulse of the thrusters. Reaction thrusters are assumed to be 
fixed along the x,y,z body axes to produce torques about the y-,z-,x-axes, 
respectively. Without loss of generality, 9 is restricted to 0 < 0 < ir/2. 

The relationships between inertial rates and body rates U)jj, Wy, and Wg, 
as determined from figure 1 (a) , are 


4 > 

0 

i 


(i)y sin (|) - a)j{ cos cf) 
sin 9 

= Uy cos (|) + o)jj sin (j) 
= Wg - 4^ cos 0 





( 1 ) 


The angular momenta along the body axes (from fig. 1(b)) are 


H cos 0 = 12^2 


(2a) 


H sin 0 sin <j) = ^y^y 


(2b) 
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-H sin 9 cos 4) = 


(2c) 


Combining equations (1) and (2) yields the Euler rate equations governing 
the coning motion, which are 










sin 0 sin cos 


* 


COS 0 

H 

Iz 


- ^ cos 0 




(3) 


A two-impulse coning maneuver is performed as follows: The body, initially 

at rest in inertial space, is acted upon by an external-control-torque impulse 
of negligible duration with predetermined components along the x-, y-, and 
z-axes. This causes the body to spin about its z-axis and precess about the 
Z-axis with nutational motion. After a given time during which the body freely 
rotates to the desired inertial attitude, a second control impulse terminates 
all motion. 

Figure 2 illustrates the maneuver geometry and shows the axial and trans- 
verse momentum components labeled ^ 2 .,o ^xy,o generated by the initial 

impulse. Also, the momentum components which terminate the motion are shown 
as H2 f and H^y f. The initial body-axis system ^o,Yor^o included to 
illustrate the orientation of H in body coordinates. 

Candidate coning solutions to equation (3) must include maneuvers with 
negative precession rates < 0) as well as maneuvers with positive precession 
rates. Also, solution paths can be classified as long or short, depending upon 
whether the total change in precession angle is more or less than tt radians. 
Four alternative solution paths (nutational motion suppressed) are illustrated 
in figure 3 for a reorientation example. 

The spherical geometry of a generalized reorientation maneuver R-j , R2r 
R3 is shown in figure 4. By definition, a R-| , R2/ R3 maneuver is any 

specified set of sequential Euler angle rotations about the z,y,z body axes, 
respectively, which reorient the body from its initial attitude to some final 
desired attitude. 

Parameter relationships for the coning motion were determined from the 
geometry of figure 4. For the spherical triangle with sides 9^, R2f and 0f, 
reference 5 specifies the following relationships: 
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Napier ' s analog i es 




sin 

1(R2 - 9o) 



sin 

"1 

2 (^2 ®o) 



”1 ^ 1 

COS 

;<R2 - Oo) 



cos 

2 (^2 ®o) 


cot 


1 

-(TT - <1)0 - Ri) 


cot 


1 

- (TT - 



(4a) 


(4b) 


Half-angle formula 


tan 




]jsin (s - 9 q) sin (s - R2) sin (s - 9f)/sin (s) 


sin (s - 9f) 


(5) 


where 


1 

s = — (9 q + R 2 + 9f) 


Law of sines 

sin 9f sin 9 q sin R2 

sin (it - ())q - R-] ) sin (<()f - R3) sin 4)f 


( 6 ) 
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Gauss ' s formula 


cos 


[^' - - 


Ri) 


cos (9f/2) sin “('J'f + 4»f “ R 3 ) 


1 

cos "{R 2 “ ®o) 


(7) 


Equations (6) and (7) were used only to check solutions. 

In addition to the geometric relationships, certain dynamic constraints 
were utilized. For example, the time between initial and final impulse (maneu- 
ver time or precession time) can be expressed by 


i|;f - 

tf - to = — : (8) 

'l^av 


where 'I'o set equal to zero by choice and is determined by 

integrating the first of equations (3) to obtain 



or-'- : ‘ ... 

• H 

- \l>av “ — (F) (9a) 


where 

Ijj + ly ly - Ijj sin (2(j>f) - sin (2(l)o) 

■ ■ 2ly 2Iy ” ^q) 


(9b) 
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Equation (9a) represents the average 4) with respect to, <}) over the coning 
maneuver and, thus, is an approximation to the average 4* with respect to time 
required in equation (8). However, the error of this approximation was found 
to be small (about 1/2 percent on the average). Combining equations (8) and (9) 
gives 


HtfF 

« 

Ix 


( 10 ) 


Another dynamic constraint to the motion, derived in reference 4, defines 
the maximum half-cone angle ^max terms of the inertias and Euler angles 
0 and 4>* For polhodes about the z-axis. 


9max = tan-1 


^y(^z tan ({) 


Ij5(Ix — ^v1 / (tan^ 0 cos^ <J)) 


(11a) 


and for polhodes about the x-axis. 


9max = tan-1 


Iy(Ix "" ^z^^x tan (J) 

|ljj(Iz — ly) ^z ^^x — ly)/ (tan^ 0 cos ^ 4^) 


(lib) 


The term "polhode" is discussed in appendix B. Equations (11) derive from the 
fact that energy and momentum of the body remain constant throughout the motion. 
These equations were used not only to insure that tentative 0,4> combinations 
produced real values for 0n,ax tut also to relate dynamically the initial and 
final 0,4> combinations through a constant ^max* 

Returning now to equation (3), it is noted that the Euler rate equations 
are nonlinear and have no known analytical solution involving only elementary 
functions. However, reference 6 formulates these equations (with the 4* and 
4> symbols interchanged) along with a solution in terms of Jacobian elliptic 
functions. The required functional relationships are 


sn (u,m) 


sin 0 sin 4> 


Q 


(12a) 


sin 0 cos 4> 

cn(u,m) = (12b) 

P 
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dn(u,m) = 


cos 0 


(12c) 


where u is the linear function of time 


u = Uq + Xt 


(13) 


and the constants P,Q,R,X, and the parameter m can be written as 


P = 


Q = 


R = 


1^- 

j - (2 TIz/h2)] 

1 Ix - Iz 

ry 

j - (2 TIz/H^)] 

11 ly - Iz 

Iz 

|]2tIx/h 2) - l] 




(14a) 


(14b) 


(14c) 


X = H 


(ly - Iz) [(2 TIx/h 2) - l] 


1-x^y^z 


(14d) 


dx - Iy)[l - (2 TIz/h2)] 

m = z = (Modulus) ^ 


(ly - Iz) [(2 TIx/h 2) - l] 


(14e) 


The energy T and angular mcmentura constants of the motion are given by 


T = -dxWx^ + Iy<^y^ + 


h 2 = dx^x)^ + dyWy)^ + dz“z) 


( 15 ) 


and the motion constant 2 T/h 2 which occurs in equations (14) can be expressed 
through equations (2) and (15) as 
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2T - fcos^ <t> sin^ A cos^ 0 

— = sin2 0 + + (16) 

h2 \ Ix ly / Iz 


It can be seen from equations (2) and (12) that the body rate solutions are 
directly related to these Jacobian elliptic functions as follows: 

(Oy 

“z 

Also from equations (12), 

0 = cos“l [r dn(u,m)] (18) 

and 




Equations (12) to (19) were not used to generate time histories of body rates 
or Euler angles, but rather served to derive Htf for equation (10). 

Equation (13) was evaluated at tf and combined with equation (14d) to 

yield 


Uf - Uq 

^ Idy - Ij) [(2 TI;,/h 2) - l] 

i X,XyX. 

Here, the quantity Uq was calculated by means of the arithmetic-geanetric mean 
process and successive use of the descending Landen transformation (both dis- 
cussed in reference 7) from sn(UQ,m) which in turn was determined from the 
first of equations (12) evaluated at tg. The value of Uf was determined sim- 
ilarly frcxn sn(u£,m) which was obtained from the first of equations (12) eval- 
uated at tf. Finally, equations (10) and (20) are combined to generate the 
dynamic expression for 4»f: 
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(Uf - Uq)F/Ix 


( 21 ) 


'I'f = 


j(Iy - Iz) [(Ix2T/h 2) - l] 




Solution of Equations 

Equations (4a), (4b), (11a) or (lib), and (5) which govern variables ipf, 
(j)f, 0f, and 4)o, respectively, were solved in a Newton-Raphson iterative 

(inner) loop in the order specified for a given value of 0 q. Starting and 
stopping conditions are given subsequently. This result was then used in a 
second Newton-Raphson iterative (outer) loop involving the dynamic value of 
(from eq. (21)) and the relationship 


®o,new “ ®o,old 


^f 

dynamic 


vO.l 


( 22 ) 


The inner loop equations were then solved again (by iteration) for the new 
value of 0Q from e<jjation (22) . This process continued with inner loop con- 
vergence for successive "new" values of 0 q until the outer loop converged and 
all equations were satisfied simultaneously. The convergence condition for the 
inner loop was that input and output values of 0 q be equal within ±0.0001 per- 
cent. For the outer loop, convergence was defined as input and output values 
of 0Q being equal within ±0.01 percent. 

The values of 0 q and iJ)q used to start the iterative process were deter- 
mined from 


‘t’o, start = IT - Hi - ^ 


(23a) 


0 


o, start 


tan“^ 


tan (R 2 / 2 ) 
cos $ 


(23b) 


Equation (23a) is an expression for the initial base angle $ of the spherical 
triangle (see fig. 4) , and equation (23b) is a geometric relationship for a 
spherical triangle initially assumed to be equilateral. 

Solutions were attempted for 16 separate values of $ (spread over the 
range 0 < 4> < tt/ 2) for each solution path. Possible solution paths (for 
each given maneuver and inertial distribution) included all eight combinations 
of long and short paths, positive and negative precession rates, and polhodes 
about the x- and 2 - axes. All such solutions were determined if possible and 
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the one associated with minimum maneuver costs was retained for computing the 
"average cost" of reorientations presented for each maneuver range. The list- 
ing and documentation of computer program IMP2, used to determine these solu- 
tions, is presented in appendix A. 


REORIENTATION COSTS 
General Considerations 

As mentioned in the "Introduction," a cost function equal to the product 
of total maneuver impulse (proportional to fuel consumed) and total maneuver 
time (nondimensionalized by spacecraft inertia about the intermediate principal 
axis) is determined for the various solution paths available for each maneuver. 
The desired or optimum path is defined as the path associated with the smallest 
cost function for each coning maneuver. Since total impulse and maneuver time 
are inversely proportional for impulsive maneuvers, results are independent of 
each of these quantities. 

It is assumed that all reorientations within a given maneuver range R 2 max 
are equally probable. Therefore, reorientation costs were sampled within each 
given maneuver range by systematically computing optimum costs for combinations 
of Rf , R 2 , and R 3 at equally spaced intervals throughout the sample space 
limits -TT ^ Rl S 7T, 0 $ R 2 ^ R2,max» -tt ^ R 3 $ it where R 2 max = 1^/360, 

V4, TT/2, 3TT/4, and n radians. These optimum-cost functions were then weighted 
(according to likelihood of occurrence) and averaged as follows: 


(Average cost) = 


^ Qoptimum cost) x sin r£[ 


2 


R2) 


In addition, the standard deviation of the optimum-cost functions was deter- 
mined with the relationship 


(Standard deviation) = 


^ Qoptimum cost) ^ x sin 


^ (sin R 2 ) 


- (Average cost) ^ 


Cost Equations 


The impulse (and thus monentum) required to initiate and terminate coning 
maneuvers has transverse exponents 


^xy,o ” ^ ®o 


^xy,f “ ^ ®f 
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and axial components 


Hz,o = H cos 9o 


^z,f ~ ^ 


as shown in figure 2. 

For reaction thrusters fixed along the x, y, and z body axes to produce 
torques about the y-, z-, and x-axes, respectively, the cost function 
Qxotal impulse) x (Total tirae)J/Iy for the two-impulse coning maneuver is 


Q2 = Q2,a + Q2,t 


(24) 


where the axial and transverse components are, respectively. 


iHltf 

Q2 A “ ®o ®f) 

lH|tf 




(25) 


Q2,T “ ‘I’ol 1^°^ ‘t’ol) 9f ^|sin (l)f ] + |cos <j)f |j 

ly L -‘J 


By ccxnbining equations (10), (24), and (25), the cost function can also be 
expressed as 


Q2 


/. . , ,\ /, . . 

= sin 0Q(|sin 4>ol + |cos (J)Q|j + sin 9f f|sin (l)f| + |cos 4>f 

lyF L 


+ cos 


9q + cos 



(26) 


If the thruster pair fixed on the spin axis to produce torques about the x-axis 
is allowed to swivel about the spin axis, the components of Q 2 t can be mini- 
mized (reduced about 21 percent on the average) at the expense of a relatively 
small thruster prepositioning term Qp with the result expressed as 

|H|tf 

Q 2 sw “ (sin 9 q + cos 9 q + sin 9f + cos 9f) + Qp 
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RESULTS AND DISCUSSION 


The results of this investigation are primarily the costs of attitude 
reorientation by the two-impulse coning method. These are presented in the 
form of average cost (over a given maneuver range) as functions of 1 - (l^/Iy) , 
Ix/Iyr and R2,max bodies of all possible inertial distribution which are 

identified in reference 4 and presented in figure 5. Although costs for the 
thousands of required individual reorientations could not be presented, such 
costs can be calculated, if desired, by means of the computer program presented 
in appendix A. 


Each of the hundreds of two-impulse maneuver solutions used to generate 
average costs for individual inertial configurations was first determined and 
then independently checked by numerical integration of the motion equations to 
assure that the desired maneuver was accomplished within an acceptable toler- 
ance. In this respect, maximum errors in ipf, proportional to R2 jnax» were 
found to be about ±0.5 percent on the average. Corresponding errors in 0f 
and (j)£ were less than ±0.005 percent. Two— impulse coning solutions for gen- 
eral reorientations of asymmetric bodies were anticipated in reference 4 on the 
basis that the three input quantities of the maneuver (coordinate-axis torques) 
should allow control of three output quantities (final Euler angles) regardless 
of inertial configuration. This is confirmed by the results of the present 
investigation. Figure 6 shows the effect of inertia distribution on average 
cost of reorientation for maneuver ranges of tt/360, tt/4, tt/2, 3h/4, and 

IT radians . The curve labeled Ijj/Iy = 1 on each plot represents the family 
of axisymmetric bodies. Most of these results were reported in reference 4 
and are included here for comparison with the asymmetric results of the pres- 
ent investigation which are the remaining curves of figure 6. 

The axisymmetric results indicate an almost linear reduction in reorienta- 
tion costs as bodies become more elongated. This is true for all maneuver 
ranges. Reorientation costs are also seen to be proportional to maneuver range 
for axisymmetric inertial configurations. 

Ins^ction of all the results of figure 6 shows that introducing inertial 
asymmetries to axisymmetric bodies causes the average cost of reorientation to 
decrease significantly for relatively "thick" bodies. (See left side of 
fig. 6.) Conversely, for relatively "slender" bodies, introducing inertial 
asymmetry results in somewhat larger reorientation costs, in retrospect, this 
result could have been anticipated because x-axis inertia is the only difference 
between comparable asymmetric and axisymmetric bodies. Since average cost is 
proportional to the effective transverse inertia of a body (a value between Ijj 
and ly) , the increase in I^ required to change the "slender" bodies from 
axisymmetric to asymmetric resulted in higher costs. Similarly, the decrease 

in Ijj required to change the "thick" bodies caused a reduction in average 
cost. 


The results of appendix C also are illustrated in figure 6. Reorientation 
costs for rod-shaped bodies with axial symmetry along the z-axis as well as 
costs for rod-shaped bodies with axial symmetry along the x-axis are in good 
agreement with both the axisymmetric and asymmetric results previously 
discussed. 
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Figure 7 was prepared to illustrate the effect of R2 max o*' average cost 
for asymmetric bodies. These curves are restricted to Ix/^Iy = 1.2 but should 
be representative of other ratios. As would be expected, average cost is seen 
to increase with maneuver range. 

The average-cost results presented in figure 6 for asymmetric bodies repre- 
sent minimum-cost reorientations involving a mixture of polhodes about the x- 
and z-axes . Although polhode axis for minimum maneuver cost is unpredictable 
for individual reorientations, the polhode axis, on the average, seems to be 
the principal axis about which inertia is most unlike (greatest percentage dif- 
ference) the y-axis (intermediate) inertia. More specifically, an inertial fac- 
tor Pjj for estimating the probability of polhode about the x-axis has been 
found empirically for R2,max = '’f radians to be 


Px = 


Rx/y 


- 1 


^x/y 


- 1 + 


Ry/z 1 


where 


^x/y “ ^x/^y ^y/^x 

^y/z “ ^y/^z ^z/^y 


(27) 


(Whichever is greater) 
(Whichever is greater ) 


It follows that the inertial factor for estimating the probability of 

polhode about the z-axis is ?z - ^x* Figure 8 illustrates the correlation 
of Px with measured values for all results of figure 6(e), 

In reference 4, reorientation by means of impulse coning was solved under 
the constraint that the initial and final angles between spin axis and the 
momentum vector be equal. This resulted in solutions involving a two-impulse 
coning maneuver followed by a spin-correction maneuver. As a result, costs 
were considerably higher than those of the present investigation which uses two- 
impulse coning maneuvers not requiring a follow-up spin maneuver. 

Maneuver costs have been presented thus far in the form of average cost 
for specified maneuver ranges. Also impor.tant is the measure of how costs are 
dispersed about these average values which is given by their standard devia- 
tion. Ratios of standard deviation to average cost are computed along with the 
average-cost results. These values are presented in table I along with compan- 
ion values of average cost. 

Table I indicates that average costs have standard deviations of about 
40 percent (of cost) for R2,max ~ radians, increasing to about 65 percent 
^2, max reduced to tt/ 360 radians. Thus, average-cost results do not 

^PPly accurately to individual reorientations, but rather to a large number of 
uniformly distributed reorientations. 
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CONCLUDING REMARKS 


This report presents data figures which can be used to evaluate the costs 
of attitude reorientation by two-impulse coning for bodies of all possible iner- 
tial distribution and over the complete range of maneuver angle. Two- impulse 
coning is a method of general-attitude reorientation in which torque impulses 
are used to initiate and later to terminate free precessional and nutational 
motion of a body. 

Costs are presented for selected inertia configurations in the form of 
average cost (and standard deviation) for a large number of uniformly distri- 
buted reorientations within a given maneuver range. These results, determined 
for five maneuver ranges, are intended primarily for cost comparisons of the 
two- impulse method with other reorientation methods (not included) . Although 
costs for individual reorientations are not included, such costs can be deter- 
mined, if necessary, through the use of an included computer program. 

Results show that the two-impulse coning method is capable of producing 
general reorientations of asymmetric as well as axisymmetric bodies. 

Results for axisymmetric bodies indicate that for all maneuver ranges 
reorientation costs become progressively smaller as the body is elongated. 

Also, reorientation costs are proportional to maneuver range. 

The introduction of inertial asymmetry to axisymmetric bodies causes 
reorientation costs to decrease significantly for relatively "thick" bodies. 

For relatively "slender" bodies, however, the addition of inertial asymmetry 
results in somewhat larger reorientation costs. In all cases, however, costs 
are proportional to maneuver range. 

Polhode axis for minimum maneuver costs is unpredictable for individual 
reorientations. However, for a maneuver range of n radians, polhode axis, on 
the average, seems to be the principal body axis about which inertia is most 
unlike the intermediate principal-axis inertia. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 25, 1979 
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APPENDIX A 


PROGRAM IMP2 

Computer program IMP2, used to determine minimum costs for two-impulse 
reorientation of asymmetric spacecraft, is presented in this appendix. 


Symbols 

Pertinent symbols used in program IMP2 and their definition or relation- 


ship to 

symbols used throughout the report are as follows; 


ACOST 

weighted average cost of maneuvers within given R 2 

^max' 


^ (Q 2 I sin R 2 I ) 



^ (| sin R 2 I ) 


ANG 

initial amplitude of ANl used to start iterative 

solution 

AN! ,AN2 

base angles of spherical triangle, opposite sides 
respectively 

0 f and 0 Q, 

B 

indicator for positive (B = 0) and negative (B = 1) 

precession rates 

EOH 

2t/h2 


EPS 

“o 


HOLD 

solution terminated for HOLD = 1 


HTF 

Htf 


LOH 

X/H 


M 

m 


PrQ/R 

P, Q, and R, respectively 


PATH 

indicator for short (PATH = 0) and long (PATH = 1) 

solution paths 

PI 

TT 


R1 ,R2,R3 

R-| , R 2 , and R 3 , respectively 


R1C,R2C, 

R3C maneuver about x polhode axis equivalent to 

maneuver about z polhode axis; R2C = R 2 c 

Rl , R 2 , and R 3 

R4,R5,R6 

ly/ and I 2 , respectively 
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R7 

R8 

RIO 

R13 

R14 

R15 

R19 

R20 

R22 

R23 

R26 

R30 

R34 

R40,R41 

R42,R43 

R44,R45 

R54 

S 

TCOST 

TNO 

UF 

UO 


®max 

4>f 

p 

9f 

latest iterative value for (})q 
maneuver cost 

latest iterative value for ANl 
AN! check value from equation (7) 
latest iterative value for 9 q 
dynamic 

lower and upper sample limit for R1 , respectively 

lower and upper sample limit for R2, respectively; R43 = 

lower and upper sample limit for R3, respectively 

sin R2f weighting factor in cost calculation 

s 

weighted sum of maneuver costs 
weighted sum of maneuvers 
Uf 


^2 , max 
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Program Listing 


PROGRAM IMP2 ( INPUT, OUTPUT ,TAPF5=IMPUT,TAPF6=nUTPUT) 

EXTERNAL EULER 

COMMON /EULER/CA,C5,C6,H 

COMMON B,PI ,R1,R2,R3,RA,R5,R6,R7,R8,R10,R11 ,R12,R13,R14,R15, 
1R19,R20,R22,R23,R30,R31 ,R32 ,R33,R34,R53,ANl , AN2 , FLG2 ,R73 , S 
1 ,R61,H0LD,K0UNT, ADD,R123,R26,RFT, AMG,HTF,R0,EnH,R34R, ACCtCTAVG 
EQUIVALENCE! IX,R4),(IY,R5),(IZ,R6) 

REAL IX,IY,IZ 
INTEGER CSAV 
INTEGER A,C,A1,CA1 
DIMENSION WK<21) ,Y(3) 

DATA A,C/1HA,1HC/ 

PI=2.«ASIN( 1 . ) $ CON=l. $ IY=100. 

NAMFLIST/ICONE/IX, lY, IZtCON, PTDEN,R40,R41 ,R42,R43,R44,R45 
PTDEN=12./PI t CON = 2. t> R40=R42 = 0, $ R44 = -PI * Ra1=R43=R45 = PI 
1000 REA0(5,IC0NE) 

IF(E0F(5) )14,15 

14 STOP 

15 CONTINUE 

R70=REP=TNO=TCOST=TC2=TCREP=R123=FLG2=0. 

ER1S=0. 

C22=1.E8 i ISAV=NUM=0 
R0=.5 $ ACC=l.E-6 $ R34=100, 

C INITIALIZE MANEUVER 

R1=PR1=S51=R40+1 ./ ( 2.«PTDEN ) 

R2=PR2=R42+l./( 2.*PTDEN) 

R3=PR3=S53=R44+C0N/ (4.#PTDEN ) 

WRITE(6,IC0NE) 

R54=ABS(SIN!R2) ) $ GO TO 2 
C NEXT MANEUVER 

7 R3 = PR3=R3+CON/PTDEN $ I F ( R3 .GT .R45 ) GO TO 25 $ GO TO 2 

8 R1=PR1=R1+1./PTDEN $ I F ( R 1 . GT . R4 1 ) GO TO 24 $ GO TO 2 

9 R2=PR2=R2+1./PTDEN $ R 54= AB S ( S I N ( R2 ) ) $ I F ( R2.GT.R43 )G0 TO 26 

C COMPUTE EULER COSINE MATRIX FOR MANEUVER R1,R2,R3 ABOUT Z Y Z AXES 

2 SR1=SIN(R1) $ SR2=SIN(R2) $ SR3=SIN(R3> 

CR1=C0S(R1) $ CR2=C0S(R2) $ CR3=C0S(R3) 

B1 1=CR1#CR2*CR3-SR1#SR3 $ B 12=SR 1*CR2#CR3+CR 1*SR3 
B13=-SR2*CR3 t B21=-CR 1=)'CR2*SR3-SR 1>!'CR3 
B22=-SR1#CR2#SR3+CR1#CR3 S B23=SR 2*SR 3 
B31=SR2#CR1 $ B32=SR2*SR1 i B33=CR2 
C COMPUTE EULER ROTATIONS RIC, R2C,R3C ABOUT -X Y -X AXFS 
R1C=ATAN2( B12,-B13 ) $ X3=-B33*S IM( R 1C )-B32*COS ( RIC ) 

Y3 = B22*C0S(R1C)+B23#SIN(R1C ) $ X2 = B2 1«X3+B31*Y3 * Y2 = B11 
R2C=ATAN2(X2,Y2) i R 3C= ATAN2 ( X3 , Y3 ) $ ANG0=.01 t ALT=0. 

3 IF( REP.EO.O. )G0 TO 1 $ R1=-R1C $ R2=R2C $ R3=-R3C 
C INITIALIZE FOR NEWTON RAPHSOM ITERATION METHOD 
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1 ANG=ANGO $ B=l. $ IFUSAV.GT.O )B = 0. 

12 H0LD=RET=0. $ K0UNT=0 $ R31=R32=R33=-10. 

R8=R61 = PI-Rl-( 2.*B-l . ) ♦ ( I .-2 . »!‘R 1 23 ) « ANG 

R7=ATAN(S0RT(RA*(R6-R5)/( ABS(Rft-R5) /(R6-R5)+Rft#(R5-R4)#C0S(R«) 
l*COS(R8)))) $ IF(R7.LT..4#R2)GO TO 18 
R31=ATAN(TAN(R2/2. )/COS( ANG) )SIF( ALT.EO.l. )R31 = (R7+R31 ) /2,.tR7 = R31 
CALL THETA i I F ( HOLD.GT. . 5 )GO TO 18 
R51=R30 S R7=R32=(R7+R30)/2. 

CALL THETA $ I F ( HOLD.GT. . 5 )GO TO 18 * R52=R30 

5 R90=(R52-R51 )/(R32-R31 ) S R7=(R51-R90*R31 ) /( 1 .-R90 ) 

R33=R7=ABS( ASIN(SIN(R7) ) ) 

CONS=l. 

21 CALL THETA S I F ( HOLD.GT. . 5 )GO TO 18 $ R53=R30 
R34A=ABS(R34R ) J, I F ( R34A .GT . P I ) R34A = 2 . *P I-R34A 
IF(ABS( (R34A-ABS(R13) )/R34A).LF.100.*ACC)G0 TO 6 

IF(ABS(R32-R52).GT.ABS(R33-R53) )GO TO 22 * I F ( COMS .GT . 1 0. ) GO TO 22 
R33=R7=R0#R33+( l.-R0)*R32 $ C0NS=C0NS+1. i GO TO 21 

22 R31=R32 $ R32=R33 $ R51=R52 * R82=R53 t GO TO 5 

6 R8=R20 $ R12=(R13+R14-R3)/2. $ R7=R30 
C SOLUTION GEOMETRY CHFCK (GAUSS FORMULA) 

R26=C0S(R19/2. ) «S I N( R 1 2 ) /COS ( .5*(R2-R7) )t IF(ARS(R26).GT.1.)R24=1. 
R26=AC0S(C0S(2.#AC0S(R26) ) ) *ABS ( R 13 ) /R 1 3 

IF{ ABS( (R23-R26)/R23).LE..001 )G0 TO 51 * WR I TE ( 6 , 4 ) R2 3 . R 26 
C COMPUTE MANEUVER COST 

51 R21 = 2.«R4*R5*R34/(R4+R5+(R5-R4)*( S IN( 2.*R14)-SIN( 2.*R20) )/( 2.*{R14 
1-R20) ) ) 

R22=ABS(R21 )*( SIN(R7 )*( ABS(SIN( R20) ) +ABS ( COS ( R20 ) ) ) +S I M ( R 1 9 ) * ( 
1ABS(SIN(R14) )+ABS(C0S(R14) ) ) +COS ( R7 ) +COS ( R 19 ) )/R5 
C POLHODE CONDITION DETERMINED 

A1=A $ IF( (REP-.5) .LT.O. )A1=C 
AN2=R14-R3 S AN1=PI-R1-R8 

19 IF(ABS(AN2).LT.PI)G0 TO 20 i AN2=AN2-2 .#P I #ABS ( AN2 ) / AN2 S GO TO 19 
C SOLUTION GEOMETRY CHECK (SIN RULE) 

20 F40=SIN(R13)/SIN(R2 ) S F41 = S I N { AN 1 ) /S I M ( R 19 ) $F42 = S I M ( AN2 ) /S I N ( R7 ) 
I F ( ABS ( { 2.*F40-F41-FA2 ) /F40) . LF. . 007 )G0 TO 52 
WRITE(6,4)F40,F41,F42 

C LOWEST COST SOLUTION IS SAVED 

52 IF(R22.GE.C22 )G0 TO 18 

CA1=A1 $ CAN1=AN1 $ CAN2=AN2 i CRFP=RFP S CCR1=R1 $ CCR2=R2 
CCR3=R3 $ CSAV=ISAV S C7=R7 S C8=RR $ C14=R14 S C15=R15 $ CHTF=HTF 
C34R=R34R $ C19=R19 $ C22=R22 $ C34=R34 $ CB=B t C23=R23 t C26=R26 
Cll=Rll $ C12=R12 $ C13=R13 S KOUNTC=KOUMT $ C123=R123 JAM=AMG 
C4=R4 i C5=R5$ C6=R6 
STAVG=CTAVG 

C CHANGE STARTING POINT AND COMPUTE ALTERNATE SOLUTION 
18 ANG = ANG+.2 S I F ( ANG .L E . P I /2 . )G0 TO 12 $ NUM=NI)M+1 
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C 

C 


C 


C 


IF(R123.GT..5 )G0 TO 10 $ R123=l. $ GO TO 1 

10 R123=0. $ ISAV=ISAV+1 $ I F ( I S AV, LT . 2 ) GO TO 1 * ISAV=0 
SETTING CONDITIONS FOR POLHODE ABOHT(A) AXIS OR RESETTING EOR 
POLHODE ABOUT (C) AXIS 

SR4=R4 $ R4=R6 $ R6=SR4 

IF(REP.GT..5)G0 TO 11 i REP=1. t GO TO 3 

11 REP=0. $ R1=PR1 $ R2=PR2 $ R3=PR3 

IF( ALT.EO.l. )G0 TO 28 

1F( ANG0.E0..11 )G0T023$ I F ( ANGO . EO. .01 ) ANGO= . 1 1 $ I F ( ALT . FO.O . ) GOTOl 
23 IF(ABS(C34).LE.2.=f=PI )G0 TO 17 S ANG0=.02 $ ALT=1. $ GO TO 1 
28 IF( ABSIC34) ,LE.2.*PI )G0 TO 17 

IF(ANG0.EQ.,17)G0 TO 16 $ I F ( ANGO. EO. . 07 ) ANGO= . 1 7 
IF(ANG0.E0..12)ANG0=.07 $ I F ( ANGO .EO. .02 ) ANGO= . 1 2 t GO TO 1 
17 IF(C22.E0.1.E8)G0 TO 7 

COMPUTE STATISTICALLY AVERAGE MANEUVER COST 
13 TN0=TN0+R54 $ TCOST= TC0ST+C224R 54 $ ACOST=TCOST/TMO 
TC2=TC2+C22*R54#C22 
TCREP=TCREP+CREP $ R70=R70+1. 

MATRIX TRANSFORMATION CHECK FOR EACH MANEUVER 

S20=SIN(C7) $ C20=C0S(C7) $ S30=SIN(C8) $ C30=C0S(C8) 

S1F=SIN(C13) $ C1F=CCS(C13) $ S2F=SIN(C19) $ C2F=C0S(C19) 

S3F=SIN(C14) $ C3F=C0S(C14) $ S1=SIN(CCR1) $ Cl=CnS(CCRl) 

S2=SIN(CCR2) $ C2=COS(CCR2) * S3=SIN(CCR3) S C3=COS(CCR3) 

M011=C20*C30 
M11=-S1*S3+C1*C2*C3 
MF11=-S1F«S3F+C1F*C2F*C3F 
M012=-C20*S30 
M12=-S1*C3-C1*C2#S3 
MF12=-S1F*C3F-C1F«C2F*S3F 
M013=S20 
M13=C1«S2 
MF13=C1F#S2F 
M021=S30 

M21=Cl*S3+Sl*C2=i'C3 
MF21=C1F*S3F+S1F*C2F#C3F 
M022=C30 

M22=C14C3-S1*C2>)'S3 
MF22=C1F«C3F-S1F=!<C2F*S3F 
M23=S1#S2 
MF23=S1F*S2F 

D11=M011*M11+M0124M21+M013*M31-MF11 
D12=M011#M12+M012*M22+M013*M32-MF12 
D13=M011*M13+M012#M23+M013*M33-MF13 
021=M021*M1 1+M022#M21-MF21 
022=M021#M12+M022«M22-MF22 
D23=M021*M13+M022*M23-MF23 


i M031=-S20*C30 
$ M31=-S2*C3 
t MF31=-S2F#C3F 
$ M032=S20*S3n 
« M32=S2*S3 
S MF32=S2F*S3F 
* M033=C20 
$ M33=C2 
S MF33=C2F 
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031=M031#M11+M032«M21+M033*M31-MF31 

D32=M031*M12+M032*M22+M033*M32-MF32 

033=M031*M13+M032*M23+M033*M33-MF33 

IF( ABS(011 )+ABS(012)+ABS(D13)+ABS(021)+A8S(n22)+ABS(D23)+ABS(D31 ) 
1+ABS<D32)+ABS(D33).GT..006) 

lWRITE(6,53)Dll»D12,D13t021»022,D23»031rn32,D33 
53 FORMAT( 1X9E11.3) 

16 CONTINUE 

WRITE(6,29)CCRl,CCR2,CCR3tC13tCHTF»C34R»C7»C19tCB»C14,CANl tCAN2, 
ICAl ,C34,C22 tPR 1»PR3,INT(C123) tlNKCB) ,1NT(R70) ,KOUNTC»AN 
29 FORMAT ( IXF8 .4» F7. 4, F8 , 4, F7. 3 » F7, 1 , F7. 3» 2F7. 4, 3F8. 4, F7. 3 » 1 XA 1 , F7. 3, 
1F6.2,F6.3»F7.3,2I2,I5»I3,F5.2) 

C SOLUTION CHECK VIA NUMERICAL INTEGRATION 

N=3 $ T=0. $ H=25.=»ABS(CHTF)/CHTF $ TENO=CHTF/H 
Y(1)=0. $ Y(2)=C7 $ Y(3)=C8 $ TOL=l.E~5 $ MTH=2 
GMIN=ABS(CHTF)/5.E11 i GMAX=20. $ G=ABS ( CHTF ) /500 . 

DO 35 I=l»900 

CALL RKF4(N,T»TENO,Y,TOL»EULFRtPOtMTH,GMIN,GMAX,G,WK,IERR) 

IF( lERR.LT.OSTOP 1234 $ IF( T.GE. TEND-1 . F-6 )GO TO 36 

35 CONTINUE 

36 ER1=100.*A8S( (Y(l )-C34)/Y( 1 ) ) $ I F ( A8S ( C34 ) . GT. 2. *P I ) FRl =0. 
ER2=100.«ABS( (Y(2)-C19)/Y(2 ) ) i ER3= 100 ,*ABS ( ( Y ( 3 )-C 14 ) /Y ( 3 ) ) 
ER1S=ER1S+ER1 S AER1=ER1S/R70 

WRITE (6 ,4) Y( 1 ) ,Y(2 ) »Y(3 ) t FR 1 » ER2 ,ER3 » AER 1 , STAVG ,C 15 
27 C22=1.E8 $ GO TO 7 

25 R3=PR3=S53 $ GO TO 8 

24 R1=PR1=S51 S SO=SQRT( TC2/TNO-ACOST*ACnST) $ SDA=SO/ACOST 
PR=TCREP/R70 

WRITE(6,4)R2»TN0,TCREP,AC0ST»SDtSDA,PR,AERl 
GO TO 9 

26 SD=SORT( TC2/TN0-AC0ST«AC0ST) $ SDA=ACnST/SD 
WRITE(6,4)R70,TN0,TCREP,AC0ST,S0,SDA 

4 F0RMAT(6X9E14.4) 

GO TO 1000 i END 


SUBROUTINE THETA 
C NEWTON-RAPHSON OUTER LOOP 

COMMON B,PI,Rl,R2tR3,R4TR5,R6»R7,R8,R10,Rll»R12,R13,R14»R15» 
1R19,R20,R22»R23,R30»R31 tR32,R33,R34»R53»ANl,AN2,FLG2,R73,S 
1 tR61tHOLOtKOUNT» A00tR123tR26,RET» ANG,HTF»RO, E0H,R34R, ACCtCTAVG 
REAL LA(10),LB(10),LC(10 ) ♦LOH,M,n»NUM,NO,NF 
K0UNT=K0UNT+1 $ I F ( KOUNT.GT. 1 5 ) HOLD= 1 . * ADD=0. 

IFIHOLO.EQ.l, IRETURN 

CALL PHI $ IF(FLG2.EQ. 1. )G0 TO 20 $ I F( HOLD. FO. 1 . IRETURN 
R71=R20 $ R62=R8=(R8+R20)/2. 
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CALL PHI $ IF(FLG2.EQ.l. )G0 TO 20 $ I F ( HOLD. FO. 1 . )RFTURM 
R72=R20 

15 R100=(R72-R71 )/(R62-R61 ) $ R63=R8=( R71-R100*R61 ) / ( 1 .-R 100 ) 

CONS=l. 

26 CALL PHI $ IF(FLG2.FQ.1.)G0 TO 20 i IFCHOLD.EO.l. )RETURN i R73=R20 
IF(ABS<R62-R72).GT.ABS(R63-R73) )GO TO 31 i I F ( CONS. GT. 10. ) GO TO 31 
R63=R8=R0#R63+( 1 .-R0)*R62 $ C0NS=C0NS+1. i GO TO 26 
31 R61=R62 i R62=R63 $ R71=R72 S R72=R73 i GO TO 15 

20 R15=(R4+R5)/(2.#R5)+(R5-R4)/(2.#R5)*(SIN( 2.*R14)-SIN(2.*R8) )/(2.# 
KR14-R8) ) $ FLG2=0. 

E0H=SIN(R19)«SIN(R19)*(C0S(R14)#C0S(R14)/R4+SIN(R14)*SIN(R14)/R5) 

1+C0S(R19)*C0S(R19)/R6 

R=S0RT(R6#(R4*E0H-1.)/(R4-R6) ) i 0=S0R T ( R 5* ( 1 .-R6*E0H ) / ( R 5-R6 ) ) 
P=SQRT(R4*( l.-R6*60H)/(R4-R6) ) 

LOH=SORT( (R5-R6)*(R4«E0H-1. )/(R4*R5#R6) ) 

M=(R4-R5)*( l.-R6«E0H)/(R5-R6)/(R4*E0H-l. ) $ SWITCH=0. 

LA(1)=1. $ LB(1)=SQRT<1.-M) $ LC( 1 )=SORT< M) i N=l. $ 1=2 
25 LA( I )=(LA( I-l )+LB(I-l ) )/2. t LB( I )=SORT( LA( I-l )#LB( I-l ) ) 

LCn ) = (LA( I-l )-LB( I-l ) )/2. $ N=N+1. 

IF(LC( I ).LT.1.E-6)G0 TO 30 S 1=1+1 $ GO TO 25 
30 N0=S1N(R8)/Q $ 00=C0S(R8)/P S NF=SIN( R 14) /O i DF=COS (R14 ) /P 
R8I=AINT(2.*R8/PI ) $ R14I=AI NT( 2.#R14/PI ) 

U0=ATAN2(N0»00) $ UOI=AINT( 2.*U0/PI ) 

IF(R8I.EQ.U0I )G0 TO 52 $ U0=U0+PI/2.«(R8I-U0I ) 

IF( AM0D(R8I,4. ).NE.AMOD(UOI ,4. ) ) UO=UO+ABS ( R8-U0 ) / ( R8-U0 ) *P I /2 . 

52 NUM=1. $ 1 = 1 $ PBN=»10 

35 QN=AINT(2.*PBN/PI )+ABS(PBN)/PBN 

QUAD=QN-ABS(PBN)/PBN $ I F( QN/2. . EO. AINT( ON/2. ) )OUAO=ON 
PB=ATAN(LB( n#TAN(PBN)/LA( I ) )+PBN+0UAD*PI/2. $ NUM=NUM+1. 

1=1+1 $ IF(NUM.EQ.N)GO TO 40 $ PBN=P8 * GO TO 35 
40 IF(SWITCH.E0.1.)GO TO 45 * SWITCH=1. $ EPS=PB/ ( 2 .*« ( N-1 . ) *L A ( I ) ) 
1=1 $ UF=ATAN2(NF,0F) i UFI=AINT(2.«UF/PI ) 

1F(R14I .EQ.UFI )G0 TO 53 t> UF=UF+P 1 /2 < R14I-0F 1 ) 

IF( AM00(R14I t4. ) . NE. AMOO( OFI » 4. ) )UF=UF+ABS ( R14-UF ) / ( R 14-UF )*PI /2. 

53 NUM=1. $ PBN=UF $ GO TO 35 

45 HTF=(PB/(2.*#(N-1. )*LA( I ) )-EPS ) /LOH#ABS ( R4-R6 ) / ( R4-R6 ) 
CTAVG=ABS(R#(UF-UO)/(LOH«HTF) ) 

R34=HTF*R15/R4 $ R34R=R34-2.*PI*INT( R34/ ( 2. «P I ) ) 

5 R53=R30=R7#ABS(R13/R34R)#*.l 

IF( ABS(R34R).GT.PI )R53=R30=R7*ABS(R13/( 2.#PI-ABS(R34R )))**.! 
IF(R30.LE.R2/6.0)H0L0=1. 

RETURN $ END 


SUBROUTINE PHI 

C NEWTON-RAPHSON INNER LOOP 
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COMMON B,PI ,R1 ,R2»R3tR4tR5t‘^6tR7tR8»R10,Rll,R12tR13tR14,R15» 
lR19,R20,R22,R23,R30,R31»R32,R33tR3A,R53»ANl ,AN2,FLG2,R73tS 
1 ,R61,H0LD,K0UNT» A00»Rl23»R26»RETt ANG»HTFtR0,E0H,R34R» ACC»CTAVG 
ADD=A0D+1, $ IFIAOD.lt. 20, )GO TO 10 
H0L0=1. $ RETURN 

10 R9=(R5*(R6-R4)*C0S(R8)«C0S.(R8)+R4*(R6-R5)*SIN(R8)#SIN(RR) )#SIN(R7) 
1*SIN(R7) 

R10=R9/(R6=*‘(R4-R5)*SIN(R7)*S1N(R7)*C0S(R8)*C0S(R8)+R4=«'(R6-R5)# 
1C0S(R7)#C0S(R7) ) $ I F(R10.GE,0, )G0 TO 20 
R47=ATAN(SQRT(R4*(R6-R5)/( ABS(R6-R5)/( R6-R5)+R6*(R5-R4)#C0S(R8) 
1*C0S(R8)))) $ RAN=ABS(R7*R8)/200. 

15 IF(RAN.LE..05)G0 TO 16 
RAN=RAN/2. $ GO TO 15 

16 R47=R47*( 1 .-RAN) 

IF(R7.EQ.R31 )R31=R47 $ I F( R7. EO. R32 ) R32=R47 
IF(R7.EQ.R33)R33=R47 $ R7=R47 $ GO TO 10 
20 R10=ATAN(SORT(R10) ) 

AN1=PI-R1-R8 

R11 = ATAN(SIN( .5«(R2-R7) )/(SIN( .5=«=(R2 + R7) ) «TAN ( . 5*AN 1 ) ) ) 
R12=ATAN(C0S( . 5* { R2-R7 ) ) / ( COS ( .5*(R2+R7) )#TAN( .5*AN1 ) ) ) 

R13=R11+R12 i PATH=1.-2.*R123 

IF(0. .GT.R13*PATH*( 2.*B-1 . ) )R 13 = R 1 3-P I *ABS ( R 1 3 ) /R 1 3 
R14=R13-2.*R11+R3 

5 IF(0..LT,R13*PATH*(R14-R8)#(R4-R6) )G0 TO 9 

R14=R14+2.#PI#PATH 

IF( (R14-R8 )#PATH.GT.2.#PI )R14=R14-4.#PI*PATH S GO TO 5 
9 R18=R54(R6-R4)#C0S(R14)#C0S(R14)+R4*(R6-R5)*SIN(R14)«SIN(R14)-R6 
1#(R4-R5)«COS(R14)*COS(R14)#TAN{R10)«TAN(R10) 
R19=ATAN(SQRT(R4#(R6-R5)*TAN(R10)*TAN(R10) /R18) ) 

S=.5#(R7+R19+R2) 

IF( (S-R2)*(S-R7)*(S-R19).GT.O. )G0 TO 8 S H0LD=1. $ RETURN 
8 R46=SIN(S-R7)«SIN( S-R2 )*SIN ( S-R19 ) /S I N ( S ) 

R23=ABS( 2.#ATAN( SQRT(R46)/SIN(S-R19) ) ) * ABS ( R 1 3 ) /R1 3 
R20=PI-R1-R23 

FLG2=0, $ IF(ABS( (R8-R20)/R8).LE.ACC)FLG2=1. 

4 F0RMAT(6X8E14.4) 

RETURN $ END 


SUBROUTINE EULER ( N ,T,Y ,0YDT ) 

COMMON /EULER/C4,C5,C6,H 
DIMENSION Y(3) ,DY0T(3) 

DYDT( 1 )=H*( SIN( Y( 3) )*SIN(Y(3) ) /C5+C0S ( Y < 3 ) )*COS( Y(3) )/C4) 
DYDT( 2)=H*( 1 ./C5-1 ./C4)#SIN( Y( 2) )«SIN(Y(3) )*C0S(Y(3) ) 
DYDT( 3)=C0S(Y( 2) )*{H/C6-DYDT( 1 ) ) 

RETURN $ END 
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APPENDIX B 


POLHODE AXES 

Poinsot's geometric solution, described in reference 8, equates the free 
rotational motion of a rigid body to the motion of the "ellipsoid of inertia" 
for that body as it rolls on the "invariable plane." The point of intersection 
between these two surfaces traces out a "polhode" curve on the ellipsoid of 
inertia. The axis surrounded by the polhode curve, or "polhode axis," is the 
principal body axis which processes about the momentum vector during a coning 
maneuver. The polhode axis, which can be an axis of minimum moment of inertia 
or an axis of maximum manent of inertia but not an axis of intermediate moment 
of inertia, is related to the initial motion and the inertial distribution of 
the body as follows; 


Polhode is about the x-axis if 


tan^ 9 cos 2 4) > 


IxUy - Iz) 

^z(^X ^ ^v) 




Polhode is about the z-axis if ) 


tan^ 9 cos^ ^ < 




(Bl) 


where 0 and are evaluated simultaneously at any time during the maneuver. 
Conditions (Bl) can be determined directly from equations (11) and the require- 
ment that Qjnax ^ real- valued. 

It should be mentioned that the derivation of equation (11a) was based upon 
0 being measured from the polhode axis. This condition is met in the present 
paper for motions with polhodes about the z-axis, but not for polhodes about the 
x-axis. Rather than rederive the motion equations with 9 measured from the 
x-axis, the two-step alternative used in reference 4 was used. That is, for 
polhode about the x-axis, the values of I^ and I^ were interchanged and the 
original Euler rotations about the z-, y-, and z-axes were replaced with a newly 
generated set of Euler rotations about the -x-, y-, and -x-axes designed to pro- 
duce the identical body reorientation. The new Euler set was determined by the 
method of reference 9 and implemented in program IMP2 of appendix A. With these 
changes, the coning geometry is transferred, in effect, from polhode about the 
z-axis to polhode about the x-axis. 
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APPENDIX C 


AVERAGE REORIENTATION COSTS FOR ROD-SHAPED BODIES 

The derivation of average reorientation costs for rod-shaped bodies is pre- 
sented in appendix C of reference 4. Costs for rod-shaped bodies with axial 
syiranetry along the x-axis (Ijj ^ 0; ly = i^) as well as rod-shaped bodies with 
axial symmetry along the z-axis {I 2 -^0; Ix = ly) are specified by the equation 


8 

(Average cost) = - R2,n>ean 

71 


(Cl) 


However, in application of this equation reference 4 uses incorrect values for 
^2, mean* Correct values of R2,mean' derived for rod-shaped bodies with axial 
symmetry along the z— axis, are expressed as a function of maneuver range by 


^2, mean 


^2, max “ ^2, max ^os R2,max 


1 


cos 


,max 


For rod-shaped bodies with axial symmetry along the x-axis, all coning maneuvers 
are "polhode about the x-axis" and R2,mean is determined by averaging R 2 c, 
the second angle of the transformed Euler angle set. The average value of R 2 c 
was determined numerically to be equal to ti/ 2 radians for all values of 
^2 ,max ♦ 

Values of R2,mean s^d average cost are listed as follows for pertinent 
values of R2,max snd for both inertial configurations; 



Selected values 

of ^2, max' radians 


0.00873 

77/4 

77/2 

377/4 

— 

7T 

Rod-shaped with axial symmetry about x-axis 

^2, mean • • • . 

TT/2 

77/2 

77/2 

77/2 

77/2 

Average cost . . 

4 

4 

4 

4 

4 

Rod-shaped with axial symmetry about z-axis 

^2, mean .... 

0.00582 

0.5181 

■i 

m 

77/2 

Average cost . . 

0.01481 

1 .319 


ml 

4 
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These values of average cost are presented in figure 6. 
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Cy = 100 


Ix 

Iz 

Cost at 

selected values of R 2 ^, 

radians 

, of - 


tt/360 

tt/4 

it/2 

3ir/4 

7T 




Av. 

cost 

SDA 

Av. 

cost 

SDA 

Av. 

cost 

SDA 

Av. 

cost 

SDA 

Av. 

cost 

SDA 


99 



4.79 

0.60 

5.57 

0.47 

5.91 

0.39 

6.05 

0.37 

195 

95 

5.01 

0.74 









180 

80 

3.72 

.71 

4.87 

.51 

5.88 

.51 

5.93 

.45 

6.02 

.43 

160 

99 



5.16 

.64 

5.73 

.48 

6.14 

.40 

6.28 

.37 

160 

95 

4.60 

.75 









160 

80 

3.54 

.68 

4.83 

.51 

5.85 

.50 

6.00 

.45 

6.08 

.42 

160 

60 

3.10 

.65 

4.28 

.48 

5.59 

.48 

6.13 

.47 

6.16 

.46 

140 

40 

1 .68 

.57 

3.16 

.38 

4.65 

.39 

5.82 

.45 

5.86 

.45 

120 

99 

4.20 

.86 

5.50 

.72 

6.26 

.54 

6.68 

.45 

6.80 

.42 

120 

95 

4.36 

.71 









120 

80 

3.12 

.60 

4.60 

.50 

5.72 

.46 

6.34 

.44 

6.48 

.43 

120 

60 

3.04 

.67 

4.42 

.54 

5.53 

.46 

6.27 

.43 

6.34 

.42 

120 

40 

1 .71 

.58 

3.12 

.39 

4.58 

.38 

5.73 

.42 

5.87 

.42 

120 

20 

.83 

.61 

2.14 

• 34 

3.53 

.36 

4.59 

.39 

4.98 

.42 

no 

10 

.43 

.65 

1 .73 

.34 

3.05 

.36 

4.07 

.39 

4.50 

.42 

100+ 

95 

4.00 

.76 









100+ 

80 

2.84 

.55 

4.69 

.56 

5.67 

.41 

6.53 

.37 

6.79 

.36 

100+ 

60 

2.81 

.68 

4.36 

.48 

5.54 

.43 

6.25 

.36 

6.48 

.34 

100+ 

40 

1 .8 

.60 

3.32 

.42 

4.67 

.38 

5.68 

.37 

5.97 

.36 

100+ 

20 

.84 

.62 

2.09 

.34 

3.33 

.33 

4.32 

.36 

4.79 

.39 

100+ 

10 

.41 

.64 

1 .65 

.32 

2.90 

.35 

3.88 

.38 

4.33 

.41 

100- 

200- 

2.73 

.58 

4.46 

.50 

5.61 

.42 

6.26 

.35 

6.53 

.34 

100- 

1 80 

2.73 

.57 

4.43 

.50 

5.64 

.43 

6.42 

.37 

6.65 

.35 

100- 

160 

2.88 

.69 

4.52 

.51 

5.70 

.43 

6.46 

.37 

6.76 

.36 

100- 

140 

2.74 

.55 

4.64 

.54 

5.70 

.43 

6.54 

.37 

6.78 

.36 

100- 

120 

2.94 

.57 

4.87 

.60 

6.09 

.47 

7.35 

.43 

7.50 

.41 

80 

180 

2.48 

.63 

3.78 

.51 

4.94 

.41 

5.95 

.42 

6.20 

.43 

80 

160 

2.62 

.72 

3.94 

.58 

5.05 

.44 

6.01 

.43 

6.25 

.43 

80 

140 

2.51 

.62 

3.88 

.52 

5.12 

.44 

6.04 

.43 

6.26 

.42 

80 

120 

2.58 

.61 

4.18 

.59 

5.41 

.49 

6.14 

.43 

6.30 

.41 

80 

110 

2.91 

.77 

4.47 

.69 

5.72 

.52 

6.23 

.43 

6.40 

.41 

80 

101 



5.12 

.69 

5.93 

.51 

6.40 

.42 

6.55 

.39 

60 

160 

2.23 

.74 

3.45 

.57 

5.21 

.55 

5.86 

.47 

5.98 

.44 

40 

140 

1 .97 

.89 

3.39 

.64 

5.32 

.55 

5.58 

.47 

5.66 

.45 

40 

120 



4.04 

.68 

5.44 

.53 

5.62 

.46 

5.68 

.44 

40 

101 



4.94 

.59 

5.46 

.48 

5.54 

.42 

5.61 

.41 

20 

120 



4.38 

.63 

4.90 

.48 

4.96 

.44 

5.01 

.43 

10 

110 



4.55 

.52 

4.55 

.43 

4.54 

.41 

4.55 

.42 








(a) Inertial rate components. 


Z 



(b) Angular momentum components. 

Figure 1 .- Orientation of body and inertial coordinate axes systems 
for Euler rotation sequence about body axes z,YrZ» 

respectively. 



z 



Figure 2.- Coning geometry showing momentum components. 
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Figure 3.- Basic coning-path options (short and long paths 
for both positive and negative coning rates) . Nutational 
motions are suppressed for clarity. 
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Figure 4.- Geometry of generalized maneuver R] , R 2 / ^nd R3 illustrating 

the spherical- triangle and spin-angle relationships. Positive angles are 
measured clockwise about axis of rotation. 




Figure 5.- Inertia constraints. All real bodies are in shaded areas 



(a) Maneuver range, tt/360 rad: 


Figure 6.- Effect of inertia distribution on 







(b) Maneuver range, tt/ 4 radians 
Figure 6.- Continued. 




(c) Maneuver range, tt/ 2 radians 
Figure 6.- Continued. 







(e) Maneuver range, it radians (all possible reorientations) 


Figure 6.- Concluded. 









Measured probability of polhode 
about x-axis 


Figure 8.- Correlation of inertial factor Pjj with measured probability 
of polhode about x-axis. R2,max ~ radians. 
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